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Effects of fermion-vacuum polarization by a singular configuration of an external static 
; vector field are considered in (2 + l)-dimensional spacetime. Expressions for the induced 
ON • vacuum charge and magnetic flux are obtained. 

^ ! 

O ■ 1 Introduction 

m ■ 

Effects of singular external fields (zero-range potentials) are comprehensively studied in 
^ ■ quantum mechanics (see |J and references therein). In this study, we consider the effect of 
t^j- ■ singular external fields on the fermion vacuum in quantum field theory. In contrast to the 
Schrodinger operator, the Dirac operator may be free of an explicit delta-function singularity; 
nonetheless, the problem of self-adjoint extension arises in both cases, albeit for different reasons 
(see, for example [[|). It is well known that a singular magnet ic-monopole external field leads 
to a 6 vacuum violating CP symmetry || ffl, ^ El . It will be shown below (see also ||) that 
a singular magnetic-string external field leads to a vacuum violating C symmetry. 

Our analysis will be based on representing the second-quantized fermion-field operator in 
the form 

tf(x,t) = ^ e- iE V B (x)a B + ^ e- tEt M*)b E - (1) 

E>0 E<0 
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Here, the symbol denotes summation over the discrete spectrum of the energy E and 

integration (with some measure) over its continuum spectrum; a E and ag (b E and are the 
fermion (antifermion) creation and annihilation operators satisfying anticommutation relations; 
and ^(x) is a solution to the equation 

#<Mx) = £^(x), (2) 

where 

" d 

H=-ia — - zV(x) + (3m (3) 
ox 

is the Dirac Hamiltonian in a static external vector field V(x). If the condition 

[ft{Hifj) - (fft^)t^] dQ = 0, (4) 



x 



where dQ is the volume element of a spatial region X, is met, the Hamiltonian H is a Hermitian 
(symmetric) operator acting in the space of functions defined on X. If, in addition, the spaces 
of the functions ip and ip coincide, the Hamiltonian if is a self-adjoint operator. 
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The integral on the left-hand side of (4) can be reduced to an integral over the surface dX 
bounding the spatial region X. As a result, relation (4) takes the form 



l I ftctip-dcr = 0, (5) 
Jax 



i 

lax 

where dcr is an oriented element of the surface dX, the normal to this element being directed 
outside the region X. 

A standard procedure leads to the expressions for the vacuum-charge density, 

= -^sgn(£)4(x)V> E (x), (6) 

and for the vacuum current, 



j( x ) = -^sgn(E)4(x)a^(x), (7) 



where 

. , r l,u > 
sgn(w) = 



-l,u< 

Performing integration in (6) over the spatial region X, we obtain the total vacuum charge 



Q (/) = f pdtt. (8) 
Jx 



As to a global quantity associated with vacuum current (7), it follows from the Maxwell equation 

]_d_ 

e 2 cht 



xB (/) (x)=j(x) (9) 



(e is the coupling constant) that a magnetic field is induced in the vacuum. This magnetic field 
is characterized by the field strength 

B (J) (x) = e 2 [ j(x) x rfx (j(x (oo) ) = B (/) (x (oo) ) = 0) (10) 
and by the total flux (in 2ir units) 

$^ = i- J B (I ydcT, (11) 

where dcr is an oriented element of the surface orthogonal to the lines of force of the magnetic 
field in the spatial region X. 

In this study, we consider a second-quantized fermion field in an external field generated by 
a source in the form of a singular magnetic string. If we direct the coordinate x 3 axis along the 
string, the strength of the string magnetic field is given by 

S 3 (x) = 27r$ (0) 5(x). (12) 

By x, we will henceforth mean a two-dimensional vector in the plane orthogonal to the string 
axis [x = (x 1 , x 2 )} ; the parameter <3>(°) is the total magnetic flux (in 2n units) of a string. It is 
natural to choose the gauge 

V 3 = 0. (13) 
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The two-dimensional vector potential V = (Vi, V 2 ) can then be defined as 

x ■ V(x) = 0, xx V(x) = $ (0) . (14) 

That the potential V(x) is indeterminate at the point x = of the plane is associated with the 
delta-function singularity in the strength B 3 (12) at this point. On the plane orthogonal to the 
string axis, the Dirac Hamiltonian has the form 

H = -ia r d r - ir- x a v {d v - *$ (0) ) + (3m, (15) 

where 

a r = a 1 cos p> + a 2 sin ip, oif = — a 1 sin ip + a 2 cos ip, (16) 
and where we introduced the polar coordinates 

r = [(x 1 ) 2 + (x 2 ) 2 ]? , ip = arctan 

This article is organized as follows. In Section 2, we determine the complete system of 
solutions to the Dirac equation in the field of a singular magnetic string. In Section 3, we 
consider the vacuum charge induced on the plane orthogonal to the string axis. In Section 4, 
we analyze the vacuum magnetic flux through this plane. The results obtained in this study 
are discussed in Section 5. Some technical details concerning the derivation of basic relations 
are described in Appendices A and B. 




2 Solving the Dirac equation in the field of a singular 
magnetic string 

It is well known that, in (2 + l)-dimensional space-time, the Clifford algebra does not have a 
faithful irreducible representation; instead, it has two nonequivalent irreducible representations. 
Accordingly, the matrices a 1 , a 2 and (3 admit the following realizations in terms of square rank- 
two matrices: 

a 1 = -a 2 , (3 = a 3 , a 2 = sa 1 , s = ±1. (17) 

Here, a 1 , a 2 and a 3 are the Pauli matrices, and the two possible values of the parameter s 
correspond to the two nonequivalent representations. 

A solution to the time-independent Dirac equation (2) with the Hamiltonian H in the form 
(15) is given by 

^b(x) = <x|^,n>, (18) 

neZ 

where 

f n (r, E)e m * 
9n(r, E)e l 

and Z is the set of integers. The radial wave functions /„ and g n satisfy the system of equations 
(-d r + r-hX)^ = (E + m)g n , [d r + r'^sX + l)\g n = (E - m)f n , (20) 



<x\E,n>=( )> ( 19 ) 
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where A = n — In the case of 7^ n' (where n' E Z), two linearly independent solutions 
to (20) that correspond to the continuous spectrum, \E\ > \m\, can be represented as 



fJ±\r,E) \ ( fW(E)J ±sX (kr) 
g n ^(r,E) J \ ±gM(E)J ±{sX+1) (kr) 



(21) 



where k = V E 2 — m 2 , and J^{z) is the Bessel function of order \i. It can be seen from (21) that 
a solution that is regular at the point r = can be chosen for all modes with the exception of 
that which corresponds to n = no, where no is determined from the condition 

- 1 < sA Q < 0. (22) 

For this mode, either the upper or the lower component of the spinor in (19) - depending on 
the choice of the plus or minus sign in (21) - proves to be irregular at r = 0, although it is 
square-integrable. It is also obvious that, for = n' (where n' E Z)), a solution that is 
regular at the point r = can be chosen for all modes. Let us introduce the quantity 

F=\ + s[{\^\}- 1 -). (23) 

Here, {| u |} stands for the fractional part of the quantity u - that is {| u |} = u — [[u]} 

(0 < {|u|} < 1), where [[u]] is the integral part of u (the closest integer to u from below or the 

integer equal to u if it is integral itself). Taking into account the relations 

no = m (0) B + sA ° = - F > (24) 
and the condition of orthonormality for states of the continuous spectrum for sgn(E) = sgn(E'), 

S(k — k') 

d 2 x < E, n|x >< x\E', n' >= o nn >, (25) 

we obtain 

M - J_ ( Vl + mE-iJ^ikr) \ r_ s ( n _ m m _ I \ , I (26) 
for the regular modes with sX > sX (l > 1); 



fn\ _ 1 ( V 7 ! + mE- 1 J v+F (kr) \ , _ ( (0) 1\ 1 



- S (n-[[$(°)]] --)--, (27) 



9n ) 2v^F V -sgn( J E) v / l - mE- l J v _ 1+F (kr) J ' V LL JJ 2/ 2' 

for the regular modes with sX < sXq(1' > 1); and 



f no \_ 1 / y/l+mE 1 [sin vJ_ F {kr) + cos i/Jp(fcr)] 



#n / 2a/tt(1 + sin2^cosF7r) \ sgn(£)Vl - m^-^sin vJi_ F (kr) - cosi/J^ 1+F (k 



r) 
(2 



1 J It should be noted that, in (2 + l)-dimensional spacetime, as well as in any spacetime having an odd 
number of dimensions, the parameter m appearing in expression (3) tor the Hamiltonian can take both positive 
and negative values. 
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for the irregular mode (sA = s\q). The parameter v is determined by the requirement that 
Hamiltonian (15) be a self-adjoint operator. We are now going to consider this issue in some 
detail. 

In the case being considered, the condition requiring that Hamiltonian be a Hermitian 
operator (5) takes the form 



2nr ^2 ^ n9n ~ Snfn) 



= 0. (29) 

r=0 



If the modes corresponding to n ^ n are subjected to the regularity condition at r = 0, 
Hamiltonian (15) as defined on the space of these modes becomes a self-adjoint operator. The 
mode with n = n cannot be subjected to the regularity condition, because we would then be 
obliged to discard solution (28), thereby spoiling the completeness of the set of solutions to 
the Dirac equation. Thus, there arises the problem of determining the boundary condition for 
the irregular mode at r = - in other words, the problem of the self-adjoint extension of a 
Hermitian operator (precisely for this mode). This problem is solved with the aid of the Weyl- 
von Neumann theory of self-adjoint operators (see, for example, [[I], [J). Since the defect index 
of the operator defined on the space of regular functions is (1,1) in the case being considered, 
the self-adjoint extension represents a set of operators that is parametrized with the aid of one 
real continuous variable (G), and the required boundary condition for the irregular mode has 
the form 

lim (\m\r) F cos (^s® + 'jjfno = -sgn(m) lim {\m\rf~ F sinfs® + -^SW (30) 
Substituting the asymptotic form of solution (28) for r — > into (30), we arrive at 



tan z/ = sgn(m£)W- — — ( -r—r ) A(F,S), (31) 





m 1 


1 E + 


m \ 


r(i- 


F) 



\m\ 



2F-1 



where 

MF, 6) = 2^^p tan + 1) , (32) 

and T(z) is the Euler gamma function. Formulas (31) and (32) establish the relation between 
the parameters v and G. The boundary condition (30) results in that the spectrum involves 
not only a continuum but also the bound state 



ac / sin Fn ( V l + m^E BS K F ( K r)e in ^ . 

™ l J - Xl 1 + (2F- l)m- x E BS V sgn(m)^T^F T E B ^K 1 _ F (Kr)e i ^+ s ^ ' ' '' ) " >) 



TT 



where k = ^Jm 2 — E BS) K^(z) is the Macdonald function of order fj,, and the bound-state 
energy E BS i\E BS \ < \m\) is determined as a real- valued root to the algebraic equation 



m + E BS ( k 



m — E B s V \m 



-A(F, 9). (34) 



'For s — 1 and m > 0, this boundary condition was obtained in fic(| . 
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It is obvious that there is no bound state for 



< A(F, 0) < oo (35) 
and that there arises a bound state for 

- oo < A(F, 0) < 0. (36) 
The bound-state energy is zero, E BS = 0, at 

A(F,&) = -1. (37) 

We also have 

sgn(E BS ) = sgn(m) , -oo < A(F, 0) < -1, (38) 

sgn(E BS ) = -sgn(m) , -1<A{F,Q)<0. (39) 

Thus, we have constructed the complete system of solutions to the Dirac equation in the 
field of a singular magnetic string. It should be noted that for the case in which s — 1, m > 
and E > 0, the results presented in this section were first obtained in [M . 



3 Induced vacuum charge 

By using the explicit form of solutions to the Dirac equation, we can find the vacuum-charge 
density averaged over all directions. We have 

1 f 2n 

P(r) = — J ^p(x), (40) 

where p(x) is determined by relation (6). 

For the contribution of the regular modes (26) and (27) to the averaged vacuum charge p(r) 
(40), we obtain 



r = -- 



1 r°° . — ' 

- / dkk J2 \E\- l J2{(E + m)[Jl F (kr) + Jl F (kr)] + 

J ;m i i 

sgn(E) 1=1 

+(E - m)[J? +1 _ F (kr) + Jl 1+F {kr)]}, (41) 



where we combined summation over I and I'. In general, the correct procedure should involve 
introducing the regularizing factor |£^| _t (t > 0) in (41) and going over to the limit t — > 0+ upon 
performing summation and integration. However, the final result remains unchanged if, instead 
of introducing a regularizing factor, we perform first summation over the sign of the energy and 
then integration with respect to k, the integral of the contribution of each modes in (26) and 
(27) with respect to k being convergent and integration with respect to k being commutative 
with summation over I. In our subsequent calculations, summation performed first over the 
sign of the energy and then over I is followed by integration with respect to k. 
Summation over the sign of energy and over / in (41) yields 

m r°° k 

Pn EG (r) = --l dk^===[Jl F (kr) - Jl{kr)\. (42) 



4W VWJ 



rn 
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With the aid of the relation 

1 2 ro ° 



du- 



VP + m 2 7T J k 2 + m 2 + u 2 ' 

we can perform integration with respect to k. Following the substitution u = y 1 q 2 — m 2 , we 
eventually obtain 

POO 

/W» = - A / dq JL— lh^qr^.^qr) - I F (qr)K F (qr)}, (43) 
J\m\ \/q 2 - m 2 

where I^z) is the modified Bessel function of order [i. The expression coincident with (43) is 
obtained by deforming the contour of integration in the complex plane (see Appendix A). 

An alternative representation of the contribution of the regular modes to the averaged 
vacuum-charge density has the form 

— 1 /»oo 2 2 

Preg(0 = "^TT / ds exp(- S 2 ~ ^)[h- F (s 2 ) - If(s% (44) 



-1 f 00 ^2^2 

(2tt)! 7o 

which can be obtained by using the relation 



1 2 



(in exp [ — u (k + m )] 



\/k 2 + m 2 \[k 

Taking into account (33)- (39), we find that the contribution of the bound state is 

sin Fit 

p BS (r) = sgn(m)sgn(A +!)[!- sgn(A)] ^ - - _ - x 

x[(m + E BS )K F (nr) + (m - £bs)#?-fM], ( 45 ) 

where the quantity A is given by (32). 

Taking into account (28) and (31), we obtain the contribution of the irregular mode in the 
form 

1 f°° k 

PirregM = -g^ J dk ^/ k 2 + m2 { Ak2Fm \ m \~ 2F l L (+) + L (-)) J -F( kr ) + 

+Ak^ 1 - F) m\m\- 2F [(m - VP + m 2 f L {+) + (m + VP + m 2 f L^] Jf_ F (kr) + 
+2[(m + Vk 2 + m 2 )L {+) + (m - VP + m 2 )L ( _)] J_ F (£x) J F (kr) + 
+2[(m - VP + m 2 )L (+) + (m + VP + m 2 )L ( _)] J x _ F (A;r) J_ 1+F (A;r) + 
+A- 1 k- 2F m- 1 \m\ 2F [(m + VWTm 2 ) 2 L {+) + (m - VP + m 2 f L(_)]J F (kr) + 

+ A-l^(l-F) m -l| m |2F [L(+) + L( _ } ] J2 i+ir(fcr)}) (46) 

where summation over the sign of energy has been performed, and 

L (±) = [Ap- 2(1 - F W|mp 2F (-m± VP~+^) +2cosF7r + A- 1 P 2F m"V| 2F (^± v / P _ f^)]"\ 

(47) 

In Appendix A, it is shown how expression (46) can be transformed by deforming the contour 
in the complex plane to arrive at the final form 

PirregM = ^ / dq Q [h-p^K^qr) - I F (qr)K F (qr)\ - / dqx 

llX J\m\ V q 2 - m 2 7T 771 J| m | 



-2(1-F). 



-2F 



"A 



\m\ 



2F 



-sgn(m)sgn(v4 + 1)[1 — sgn(A)] 



+ 2 + A- 1 
sin Fir 



2{1-F) 



(2tt) 2 m + E BS (2F 
x [(m + E BS )K 2 F {^r) + (m - F^i^M]. 



(48 



Summing (43), (45), and (48), we find that the averaged vacuum-charge density is given by 



p(r) 



sin Fn 
7r 3 m 



dq 



\m\ 



x 



t 1 + A {wy {1 ' F) ] K Uir) - [1 + A- 1 (^)" 2 Vi 2 -f(^) 



.4 



2F 



+ 2 + A 



-i 



g 

M 



2(1-F) 



(49) 



This expression tends to zero in proportion to \m\ 2r~2 exp(— 2|ra|r) for r — > oo and diverges 
in proportion to |m|r _1 for r —>■ 0. In the case of G 7^ |(mod-7r), integration of (49) over the 
entire plane yields the expression for the vacuum charge induced by a singular magnetic string. 
The resulting expression has the form 



Q 



in 



sgn(mj 
2^~ 



dv F[l + Av~ 1+F ] - (1 - F) [1 + A" V 



y/v=T 



Av F + 2 + A- V" F 



(50) 



It should be emphasized once again that relations (49) and (50) hold only in the case of 
nonintegral values of the string flux (0 < F < 1). For integral values of the string flux (F = 0), 
the density and the flux vanish because all square-integrable modes are then regular for r — > 0, 
and the latter case does not differ in the least from the string- free case ($'°) = 0). 

Expression (50) can be reduced to the form 



Q 



U) 



■^sgn(m) 



F H — arctan 

7T 



(l + Av F ^ 


i>=oo- 


W«-i/ 





(51) 



whence we eventually obtain 

— |sgn(m)[F — sgn(A + 1)], 



Q 



< F < i 



— - ssgn(m) arctan[tan(®)], F = ~ 

k |sgn(m)[l — F — sgn(A~ 1 + 1)], § < F < 1 



(52) 



For the case of F = |, s = 1, and m > 0, the last relation was obtained in |TT|. We also have 



lim Q* 7 - 1 = -sgn(m), 
f^o 2 

lim = — sgn(m). 
f^i^ 2 y 1 



(53) 
(54) 



In deriving the last two relations, we also considered that 

lirn A = 0, ^ s|(mod27r), 

limyT^O, ^ -s^(mod27r). (55) 

F — >1 Z 

From relations (53) and (54), it can be seen that the vacuum charge as a function of the string 
flux undergoes discontinuities at integral values of its argument. 

This is confirmed in the case of = |(mod7r) as well. From (49), we can easily obtain 



sin Fir „_ f 00 



7T 3 7X1 J\m\ ^J q 2_ m 2- 

p(r) = { . (56) 

S ^ m G dq V^ K '-^ © = -*!(niod27r) 

Integrating this relation over the entire plane, we arrive at 

|sgn(m)F, = s|(mod27r) 



SZ^-rhs^M, = 4(mod27r) 



q(I) = { -|sgn(m)F, = sf(mod27r) 

" I |sgn(m)(l — F), = -sf (mod 2^)' 1 J 



4 Induced vacuum magnetic flux 

By using the explicit form of solutions to the Dirac equation, we can find the vacuum current 
averaged over all directions. We have 

1 f 2w 

j(0 = ^ J ^j(x), (58) 

where j(x) is given by (7). Having verified that the radial component of the vacuum current 
vanishes, j r (r) = 0, we proceed to considering the angular component J v (r). 
For the contribution of the regular modes (26) and (27) to j v (r), we obtain 

Perfornung summation over I, we arrive at 

= / dk ^=— [Jl- F {kr) - J 2 . F (kr)J- F (kr) + J 1+F (kr) J-i +J? (fcr) - J|(fcr)]. 

(60) 

In the same way as that used to evaluate the contribution of the regular modes to the averaged 
vacuum-charge density, the integral in (60) can be reduced to an integral featuring modified 
Bessel functions in the corresponding integral. The result has the form 

poo 2 

CgW = "A / dq J. {h- F {qr)K F {qr) - I F {qr)K x . F {qr) + 

Z7r J\m\ \/q 2 - m 2 



+- sin Fn[qrKl_ F (qr) - qrK 2 F (qr) + 2{F - \)K F (qr)K^ F (qr)}} . (61) 

71 Z 
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With the aid of (33)- (39), we find that the contribution of the bound state is given by 



jtSr) = ssgn(m)sgn(A + 1)[1 - sgn(A)] 



sin Fn 



K 



2n 2 m + E BS (2F-l) 



K F {Kr)K^ F {Kr). (62) 



Taking into account (28) and (31), we obtain the contribution of the irregular mode in the 
form 



(r) 



dk 



4tt J Vk^T 



- : {Ak 



-2(1-F) 



m\m\ 2F [{m — Vk 2 + m 2 )L( + ) + 



+ (m + Vk 2 + m 2 )L ( _ ) }J^ F (kr)J l ^ F (kr) + [L {+) + L ( _)][J_ F (Ax) J^ 1+F {kr) - J F (kr) J^ F (kr 
+A- 1 k- 2F m- 1 \m\ 2F [(m + Vk 2 + m 2 )L {+) + (m - Vk 2 + m 2 )L { _)] J F (kr) J_ 1+F (A;r)}, 

where the quantities L( + ) and L(_) are given by (47). In Appendix B, it is shown how expression 
(63) can be transformed by deforming the contour of integration in the complex plane to arrive 
at 



' IRREG 



2tt 2 

s sin Fn 1 x 



dq 



q 



\m\ 



-.[h- F (qr)K F (qr) - i" F (<?r)Ki_ F (gr)] 
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A(^) 2F - A~\^ 

v |m| ' v \m\ 



If q \2(1-F) 



(Tz F y^(^)^i-F(gr) 



sin Fix 

- SS gn(m)sgn(A + 1)[1 - sgn(A)]-^- ________ K F ( Kr )K 1 _ F (/ t r). (64) 

Summing (61). (62). and (64). we find that the averaged vacuum current has the form 



fir) 



s s in Ftt f°° 

7T 3 J\rn\ 



dq 



(f 



-.{qr[Kl_ F {qr) - K 2 F {qr)} + 



A (\m\) A 1 (\m\) 

A(^r + 2 + A-^) 



+ [2(F--) + 



2' A(JL 



l(_q_\2(l~F) 

—]K F (qr)K^ F (qr)}. 



-U 1 \*0--F)- 



(65) 



This expression tends to zero in proportion to |m| 1//2 r -3 / 2 exp(— 2|m|r) for r — > oo and diverges 
in proportion to r~ 2 for r — > 0. 

Averaging relation (10), where any point that is infinitely remote from the string (that is, 
any point lying on the circle r = oo) can be taken for x^), we find that the averaged field 
strength is given by 



B 3 {I )(r) 



s e 2 sin Fn 



h / dr' / dg J— = {g/[K 2 _ F ( g rO - K 2 ( g r')] + 

n Jr J\m\ \/q 2 - m 2 



^A) 2F + 2 + ^(A) 



+[2(F--) + 



2 7 4r 9 



1( _q_^2(l-F) 

lf q ^^KFiqr'^Air')}. 



(66) 



For the total flux (11) of the vacuum magnetic field induced by a singular magnetic string, we 
obtain 



$(-0 



se 2 F(l-F) 
2it\m\ 



i i i r 00 dv 



Av F - A~ l v 



-ln,l-F 



Av F + 2 + A^v 1 ^ 



(67) 
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The coupling constant e has dimensions of y \m\. Expression (67) can also be recast into the 
form [compare with (50)] 

w = se 2 F(l-F) /■« dv (1 + F)(l + Av F ) - (2 - F)(l + 

127r 2 |m| ^ ^v^T Av F + 2 + A~ 1 v 1 ~ F ' 1 J 

In the case of half-integer values of the string flux, we have 

e 2 6 1 

$ (/) = --— arctan[tan(— )], F — -. (69) 

It should be emphasized that, in contrast to the vacuum charge, the vacuum magnetic flux is 
continuous at integral values of the string flux. 



5 Discussion of results 

It has been shown that, on the plane orthogonal to the singular magnetic string specified by 
(12), the fermion vacuum is characterized by the quantum numbers (52), (57), and (67), which 
depend on the parameter G of self-adjoint extension. In relation to the vacuum for the case 
of a monopole || |, |[ 0], the 6 vacuum for the case of a string possesses a richer structure 
[dependence on s and sgn(m))]. In all probability, this is due to a nontrivial topology of 
the base space in the latter case: 7i"i = in the case of a space with a punctured point, and 
7i"i = Z in the case with a removed line (or in the case of a plane with a punctured point), 
where tci is the first homotopic group. It should be noted that the vacuum charge changes sign 
and the vacuum magnetic flux remains unchanged under either of the substitutions s — > — s 
and m — > —m. 

Let us compare expressions (52) and (57) obtained here for the vacuum charge with the 
expression for the vacuum charge induced by a regular configuration of an external magnetic 
field. In the latter case, one has |T^| 

1 If 

Q {1) = — ssgn(m)$, $ = - ^B 3 (x), (70) 

2 2n J 

where B 3 (x) is a function that is continuous everywhere with the exception of integrable singu- 
larities at isolated points or on isolated lines. In contrast to (70), expressions (52) and (57) are 
periodic in the flux of an external magnetic field. This can be considered as a manifestation of 
the Aharonov-Bohm effect |13| in quantum field theory (see [0). Since it is sometimes stated 



that the vacuum charge is not periodic in (see, for example, [[15], 0), we will dwell on this 
point at greater length. 

Under charge conjugation, 

C:V^-V, ij^a^*, (71) 

the charge operator and its vacuum expectation value, as well as the vacuum magnetic flux, 
must change sign, but this is not the case for expressions (52), (57), and (67), because the 
boundary condition (30) violates charge- conjugation symmetry. However, for a specific choice 
of the parameter B, this symmetry can be conserved. 
In particular, the choice of 



6 = sf (mod27r), a$(°> > 
6 = -s| (mod27r), s<$>^ < 



| ($(°)^n, neZ), (72) 
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which corresponds to the boundary condition considered in [17], |IBfl , leads to the expression 

(73) 



(compare with the results presented in |y| |16i ) 

-|ssgn(m){|$(°)|}, $(°)>0 
issgn(m)(l-{|$(°)|}), $(°)<0 



Q 



in 



{\ $ (0) 1} 1 o. 



This result changes sign under change conjugation, but it is not periodic in 

The parameter G can be chosen in such a way as to conserve both periodicity in <3>(°) and 
the discrete symmetry (71). This can be achieved, for example, by setting 



6 = sf(mod2^), -I < s ({|$(o)|} - I) < 

= O(mod27r), {$ (0) |}=§ 

6 = -sf(mod27r), 0<s({|$W|}-i) < \ 



(74) 



which corresponds to the condition of minimal irregularity, i.e., to a radial wave function that 
diverges for r — > no faster than r~ p , where p < |. It is with this boundary condition that the 
result reported in |I4[ is recovered in the form 



Q 



-s sgn(mj 



^gn ({|$( )|}-i)-{|$(°)|} + i 



(75) 



where 



sgn(-u) 
0, 



u ^ 
u = 



It should be noted that expression (75) is continuous for integral values of the string flux and 
displays discontinuities at half-integer values. 

Another choice of that is also compatible both with periodicity in and with symmetry 
(71) is 

6 = 0(mod27r), < {| $ (0) |} < 1. 



We then arrive at the expression 

1 



Q {1) 



--ssgn{m) 



|s gno ({| <^(°) 1} - |) + {| ^ (0) 1} - | 



(76) 



(77) 



which is discontinuous both at integral and half-integer values of the string flux. 

When the boundary condition (72) is used, expression (67) for the vacuum magnetic flux 
takes the form (see also f0 



12ir\m\ 



127r|m| 



|$(0)| }( 1_ { |$(0)| } 2 ); 

;i-{|$(°)|})[i-(i-{|$(°)|}r 



$(o) > o 
$(0) < o ' 



(78) 



We also have 



$(o) |}( i_{|$(0) |} 



127r|m 

for the boundary condition (74) and 
se 2 F(l-F) 



W{|* (0) |} 



1. 
2' 



{I * (0) I) 



(79) 



2n\m\ 



-(F ) + — 

6 V 2 1 An 



dv 



C F v F - C^v l ~ F 



i v C F v F + 2 + Cp l v 



JO) 
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for the boundary condition (76). 

In the last expression, the quantity F is given by (23), and 

C F = 2 1 " 2F£ ^p- (81) 

Expressions (79) and (80) are odd under charge conjugation and are periodic in the string 
flux. 

In conclusion, we note that the general form of the boundary condition that conserves both 
C symmetry and periodicity in the string flux is given by 

9 = e c (mod27r), -i< s ({$(o)|}_I)<o 

= O(mod27r), {| $(°) |} = I , (82) 

= -6 c (mod2^), < s{{\ $<°> |} - ±) < \ 

where — ir < 9c < tt. 



Acknowledgements 

I am grateful to S.A. Yushchenko for the contribution that he made at initial stages of this 
study and to L.D. Faddeev, A.V. Mishchenko, V.V. Skalozub, and V.I. Skrypnik for discussion 
on results presented here. 

This work was supported in part by the Ministry for Science and Technologies of Ukraine 
and by the American Physical Society. 



Appendix A 

On the basis of relations [see, for example, |19| 



J^iz) = exp(-fXTT)I^z), -it < arg z < -, 



7T 

•2' "' " x " " " " " ~ 2 
In(-z) = exp(iyU7r)^(^), K^-z) = exp(-i//7r)if /1 (z) - inl^z), -tt < arg z < 0, 



we can easily obtain 

1 % t 

JJkr)JJkr) = — {exp[-(/i - v)7r)lJ-ikr)KJ-ikr) - exp[-(v - fj,)ir)lJikr)KJikr) + 
lin I A 

+ exp[-(i/ - [i)T(}I v (-ikr)K^(-ikr) - exp[-(^x - v^I^ik^K^ikr)}. (A.l) 
With the aid of (A.l), expression (46) can be recast into the form 

p IRREG (r) = / du>F(u). (A.2) 

Here, u = k 2 is the new variable of integration; the contour circumvents the real positive 
semiaxis of the variable u, going along it at infinitely small distances from below and above; 
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and the integrand has the form 

a -y 

= ^-{Au F m\m\~ [L( + ) + L(_)]i"_,p(r\/— uj)K F (r^/— uS) + 

(47r) £ 

+A^~ 1+F m|m|~ 2F [(m — e) 2 L( + ) + (m + e) 2 L(_)]/ 1 _ i 7'(rv / — ^)^i-f(^V— ^) + 
+cj F [(m + + (m - £)L ( _ ) ](-o;)~' F / F (rv / ^)-K'F(^v / -^) + 

+w~ F [(m + + (m - s)L(-)}{-lo) F I_ F (r^/^u)K F (r^/^uj) + 

W~ F [(m - e)L w + (m + £)L ( _ ) ](-^) _1+F /i_ F (r v /= ^)^i-F(rv / ^) + 
+c<j" 1+F [(m - e)L (+) + (m + e)L ( _ ) ](-a;) 1_F /_ 1+F (rv /z U)Ki-F(rv /Z ^) + 
+,4~ m _ |m| [(ra + e) + (m — e) L(_)]Jir(rv— uj)K F {r\J— uj) + 
+A~ 1 o; 1 - F m- 1 |m| 2F [^(+) + ^(-)]/-i + F(rv /Z ^)^i-F(r V /Z ^)}, (A.3) 

where e = a/w + By continuously deforming the contour of integration in the complex u; 
plane as is shown in the Figure, we arrive at the relation 

/ dujF{uj) = I dujF{uj) + I dwF{u)+ I duF(uj) + / duT{uj). (A.4) 
JC(± JCxj Jc^d Jcrx Jcq 

The integrals along the semicircles C^j and of infinite radii vanish, whereas the integral 
along the contour circumventing the cut for uj < — m 2 can be represented as 

/ du,F(<*) = -TtV r "T==f Uu F m\m\-* F {e^[R§ + i£>] + 
JctD (4tt) J m a Vu-m 2 V w ^ > 

+e- lF7r [R[-] + i?f: ) ) ]}/_ F (r v ^)X F (r v ^) - A U - l+F m\m\~ 2F ' {e iF *[{m - iy/u - m 2 f ^j+j + 

+ (m + iVw - m 2 ) 2 i? ( ( +j] + e- iFw [(m - Wu - m 2 ) 2 i? { ( ;j + (m + iVw - m 2 ) 2 /?^]} x 

x/ 1 _ F (r v ^)K 1 _ F (r v ^) + {e iF "[(m + Wu - m 2 )i? ( ( +j + (m - iV« - m 2 )i? ( ( +j] + 

+e _iF7r [(m + Wu - m 2 )fl[~j + (m - iV« - m 2 ) J Rj:j]}/ F (r v ^)K F (r v ^) + 

+{e- iF "[(m + iVu - m 2 )/?^] + (m - zVw - m 2 )/?^] + e lFw [{m + i^u - m 2 )i? ( ( ;j + 

+(m - iVu - m 2 )R\tyl_ F (ry/^K F (ry/u) - {e^ F "[(m - iVu - m 2 )/^ + 

+ (m + iVu - m 2 )i?[+j] + e iJ?,r [(m - Wu - m 2 )/^ + (m + iVw - m 2 )i?[:j]} x 

x/ 1 _ F (r v ^)K 1 _ F (r v ^) - {e iFn [(m - Wu - m 2 )#[+j + (m + iV« - m 2 R { + ) ) ] + 

+e^ Fw [(m - Wu - m 2 )i?[;j + (m + iVw - m 2 ) J R ( ( :j]}/_ 1+F (r v ^)ir 1 _ F (r v ^) + 

+A- 1 u- F m- 1 \m\ 2F {e- iF7T [(m + iy/u - m 2 f R^ + (m - iy/u - m 2 ) 2 i?{+j] + 

+e iF7T [(m + Wu - m 2 ) 2 ^;} + (m - iVw - m 2 ) 2 i?[:}]}/ F (rv/^)K F (rv^) - 

-A-^^V-Vl^le" 1 ^^) + ^(-)] + e ^K;) + ^:]]}/-i+F(rv^)^i-F(rv^))(A.5) 
where 

= [AM~ 1+F m|m|~ 2F e iFw (m =p iVw - m 2 ) + 2 cos Fvr + 
+A- 1 n" F m- 1 |m| 2F e" iF "(m ± iVw-m 2 )]- 1 , 
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i?y = [Au- 1+F m\m\- 2F e- iF7T (m =f iV« - m 2 ) + 2 cos Fvr + 

+A- 1 W - F m- 1 |m| 2F e iF "(m ± ? v / «-m 2 )]- 1 . (A.6) 
Expression (A. 5) can be reduced to the form: 

/ duF{u) = -— r - = (- am Fn{Au F m\m\-™[e iF *(R™ + + 

Jet, 87r Jm 2 VM - m 2 Vl 

+e- iF7r (i?[;] + + (m + 2v / ^^)(e-^ J R{+J + e^i?{;J) + 

+ (m - * v / M-m 2 )(e^ J R ( ( +j + e lF "i? { ( :j)}A^(rv^) + 

+2m[I F (ry/u)K F (r^/u) - h^ry/^K^p^ry/u)] - -sinF7r{(m - W u - m 2 ){e iF * R { ?X + 

+e- <Fw i2j- ) ) + (m + ^ v V^ 2 )(e^i?j!j + e^i^l]) + A~ V^m" Vr[e~* F "(#(+) + 

+i?«) + e^(i?[;] + R{z])]}Kl F (rV^)). (A.7) 

As the result of further simplifications, we will arrive at the terms in relation (48) that are 
represented as integrals with respect to the variable q = \fu. 

It remained to consider the integral along the contour circumventing the pole of the function 
T{d). We have 

duT{u) = 2iri Res ^(w), (A.8) 



L 



u)=—k- 

where k 2 = m 2 — E BS 2 , and the quantity E BS is given by (34). Choosing the branch for 
fractional exponents according to the prescription 

(-k 2 Y = re 2 " exp(^Tr), < \i < 1, (A.9) 

we arrive at 

/ duF{u) = - * . [A K 2F m\m\- 2F e iFir I- F (Kr)K F (Kr) - 

JCq OTY\E BS \ 

-AK-^-^mlml-^e^imTlEBslfh-FMKiMKr) + 
+e lF7T (m ± \E BS \)I F {Kr)K F {Kr) + e~ lF *{m ± \E bs \)F f {kt)K f {kt) - 
- e - iF *{m =F |^ 5 |)/i-F(«r)X 1 _ F ( K r) - e* F7r (m =F |£ BS |)/_ 1+F (rer)AViKrer) + 
+A- 1 K- 2F m- 1 \m\ 2F e- iFn (m±\E BS \) 2 I F (Kr)K F (Kr) - 
-A- 1 /€ 2(1 - F) m- 1 |m| 2F e- iF7r /-i + F(^)K 1 _ F («r)] Res L (±) = 

' [(m ± |£ SS |)K 2 (kt) + (m =F l^ssO^M] Res L (±) , £ B5 ^ 0(A.10) 



2n 2 \E BS 
Taking into account the relation 



Res L (±) = , — r , E BS > , (A.ll) 

W =» K 2 (±) isinFyr \E BS \(2F- 1) ±m < ' V ; 

we obtain 

jf du; F{u>) = sgn{E BS f U 2 ^ m + E ^ {2 F-l) [( ™ + Ebs) ^ M + ( ™ ~ " 

(A.12) 
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Naturally, the same expression is obtained if the branch for fractional exponents is chosen 
alternatively as 

{-k 2 Y = it 2 " exp(-ifj,n) , < ^ < 1- (A.13) 

With the aid of relations (35)- (39), expression (A. 12) is reduced to the form coincident with 
that of the last term in (48). 

It should be also noted that the above method can be used to reduce expression (42) 
to the form (43), in which case the integrand naturally does not have poles on the segment 
-m 2 < k 2 < 0. 



Appendix B 

The contribution of the irregular mode to the averaged vacuum current can be found by a 
method that is similar to that used to calculate the contribution of this mode to the averaged 
density of the vacuum charge. In the case being considered, the function J-{uj) has the form 



1 



Au F m\m\ 2F [(m — £)£(+) + (m + £)£(-)] (— 



CO) 2 x 



x [Jx_ F (ry/^uP) K F (ry/^u) - I- F (ry/^u) K^ F (ry/^u)] - [L {+) + L ( _)] x 
x {cu 1 ^ (-w)~* +J? [h- F (rv^7) Kp (rv^Z7) - I_ F (ry/^u7) K ± _ F (rV=w)] + 

+cu F {-oj)^~ F [If (rV^P) #i-f (ry/^oJ) - I_ 1+F (ry/^U) K F (ry/^uV)] } + 
+A~ 1 uj 1 ~ F m~ 1 \m\ 2F [(m + e)L {+) + (m- e)L(_)] [Jj? (rv 73 ^) #i-f (rv 737 ^) - 

-J_ 1+F (rv^7) K F {ry/=^)]) • (A.14) 

The integral along the contour circumventing the cut for u; < — m 2 (see Figure) can be reduced 
to the form 



+ ( m + iy/u — m? ) 



dull-?- 

u 



a — 1 + F I I — 2F 



|e iF " (m-iy/u-m 2 ^ + 



+ e" 



-iFir 



m — iy/u — m 2 j + + iy/u — m 2 ^ 



—iFtt 



, d(+) 

(+) + ^(-) 



+ e' 



iFtt 



jF-7r 



d(+) , d(+) 



x (ry^) X F (rv 7 ^) - I_f (rv^) #i-f (rv^)] + |e 

} [h- F [ry/u) K F [ry/u) - I-f {rV^) Ki-f [rV^j] - [e 
+e~ iF * [R{-] + R\Z]] } [I F (ry/u) K X _ F (ry/u) - J_ 1+F (r^) K F (r^)] 
^^m^lmF {e~ iF7V (m + iy/u - m 2 ) + (m - zVw - m 2 ) # 

+ (m - iVw - m 2 ) | [J F (r>/u) #i-f (r>/w) - 



+e 



jF7T 



+ 



L (-) 



+ 



rn 



+ iy/u — 



-I 



l+F 



(A.15) 



where the quantity i?S is given by (A. 6). Simplifying the last expression, we arrive at the 
terms in relation (64) that are represented as integrals. 
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For the integral along the contour circumventing the pole of the function J-(oS) (see Figure) 
we have 



/ du = ^ { [AK-^ 1 - F) m\m\- 2F e iFlt (m t \E bs \) + e~ iFn ] [h„ F («r) K F 

Cq 



nr) 



KT 



-I-f («r) K^p (kx)] - [e iFw + A-^-^m-Vl^e"^ (m ± |£ s5 |)] [/*■ («r) K X _ F 
-I- 1+F («r) («r)]} Res o L (±) , £ B5 ^ , (A.16) 



CJ — — K 

where we choose the branch (— K 2 ) fl according to (A. 9). Taking into account (A. 11), we obtain 

sin Fir k 3 

duj'J- (u) = ssgn^Bs] 



f sin Ftt 

/ duT (u) = ssga(E BS ) KF^K^F^r) . (A.17) 



With the aid of (35)- (39), expression (B.4) is reduced to the same form as that of the last term 
in (64). 
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